Scattering theory of chiral Majorana fermion interferometry 
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Using scattering tiieory, we investigate interferometers composed of ciiiral Majorana fermion 
modes coupled to normal metal leads. We advance an approach in which also the basis states 
in the normal leads are written in terms of Majorana modes. Thus each pair of electron-hole states 
is associated with a pair of Majorana modes. Only one lead Majorana mode couples to the intrinsic 
Majorana mode whereas its partner is completely reflected. Similarly the remaining Majorana modes 
are completely reflected but in general mix pair-wise. We demonstrate that the charge current can 
also be expressed in terms of interference between pairs of Majorana modes. These two basic facts 
permit a treatment and understanding of current and noise signatures of chiral Majorana fermion 
interferometry in an especially elegant way. As a particular example of applications, in Fabry- 
Perot-type interferometers where chiral Majorana modes form loops, resonances (anti-resonances) 
from such loops always lead to peaked (suppressed) Andreev differential conductances, and nega- 
tive (positive) cross-correlations that originate purely from two-Majorana-fermion exchange. These 
investigations are intimately related to current and noise signatures of Majorana bound states. 



I. INTRODUCTION 

Majorana fermions (MFs) are their own anti-particles. 
They have been the subject of theories for more than 
seven decades without an experimental signature. How- 
ever recently it has been suggested that these exotic par- 
ticles can be discovered in condensed matters systems, 
in particular in topological superconductors (TSCs)^"— . 
TSCs are superconducting systems that admit excep- 
tional boundary states -owing to nontrivial topology as- 
sociated with bulk spectra- inside the quasi-particle ex- 
citation gaps. These exceptional boundary states are co- 
herent superpositions of electrons (particles) and holes 
(anti-particles) which, upon proper compositions, repre- 
sent realizations of Majorana states. The presence of 
Majorana states in TSCs is robust because of their deep 
roots in the global properties of bulk states. This leads 
to renewed efforts to reveal MFs in laboratories. 

Two examples of TSCs are chiral p-wave superconduc- 
tors in one dimension (ID) and two dimensions (2D). In 
ID, a chiral p-wave superconductor may accommodate a 
pair of Majorana bound states (MBSs) separately at its 
two ends-. The energy of this pair of MBSs at large sep- 
aration is exponentially close to the chemical potential - 
by definition the zero energy- of the superconductor, and 
each MBS hence provides a container for a MF to stay. 
Recently, systems equivalent to ID p-wave superconduc- 
tors have been proposed based on proximity effects in 2D 
topological insulatorsfiii, or even more close to experi- 
mental reality, in ID semiconducting quantum wires^^— . 
To detect MBSs in these system s ^"'^d^ , and possibly to 
manipulate the MFs therein, is of great interest in current 
researcbi^^. 2D chiral p-wave superconductors^,, and 
equivalent systems based on the proximity effect in topo- 
logical insulatorsiS or semiconductors* , are hosts for chi- 
ral Majorana modes (xMMs), which are gapless, charge- 
neutral edge excitations. xMMs can serve as coherent 
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FIG. 1: Schematic picture of a typical chiral Majorana 
fermion interferometer composed of a chiral Majorana mode 
(xMM) coupled to one or several normal metal leads. The chi- 
ral Majorana mode lives at the edge of the underlying topolog- 
ical superconductor (TSC), which has a full pairing gap and 
is grounded. The interior of the topological superconductor 
may accommodate a number of vortices, each permitting a 
magnetic flux quantum to penetrate through and containing 
a Majorana bound state. Majorana fermions traveling in the 
chiral Majorana mode pick up a phase that encodes (the par- 
ity of) the number of vortices, then scatter into normal leads. 
To generate charge current and noise Majorana fermions in- 
terfere pairwise. 



transmission channels for MFs, and hence are ideal in- 
gredients for building MF interferometersii"— (see, e.g.. 
Fig. [1]). Unlike conventional (electronic) interferome- 
ters, these MF interferometers do not permit an arbitrary 
magnetic flux insertion due to the presence of the un- 
derlying superconductors. Instead, magnetic flux in the 
MF interferometers is trapped in superconducting vor- 
tex cores and remains quantized in integer multiples of 
a flux quantum. In addition, each vortex here allows for 
one MBS^i^i^, and the parity of the number of vortices 
gives rise to a Z2 -type of MF interferometrj*ii"— . 

Both examples mentioned above belong to the same 
universality class, labeled D in the Altland-Zirnbauer 
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classification^^. This means that in these systems time- 
reversal symmetry is broken but particle-hole symmetry 
(PHS) -as far as the Bogoliubov-de Gennes theory is 
concerned- is present. Effectively these systems are spin- 
less, and the PHS operator can be simply expressed as 
S — cTxlC with K, the complex conjugate operator and Ux 
the Pauli matrix reversing the electron and hole compo- 
nents. We will limit ourselves to this universality class 
in this paper, and use the inherent PHS as the only con- 
straint to the scattering theory that will be developed 
here. 

Scattering theory in mesoscopic physics^ir— has been 
proven a very useful tool in dealing with electronic co- 
herent transport phenomena. Recent developments have 
also shown powerful application of scattering theory to 
the topological classification of condensed matter—"—, 
extending an approach known from adiabatic pumping 
processes^. In this paper we employ scattering theory 
to investigate interferometers built on xMMs. We pursue 
especially a general understanding of how MFs -being 
exotic (quasi-) particles which are charge-neutral- can 
be involved in charge transport. For such a purpose we 
adopt consistently the Majorana basis in our treatment, 
where the scattering theory of (charge) current and noise 
takes interestingly an off-diagonal form. Namely, it is the 
pairwise interference between transmitted MFs, instead 
of the self-interference of individual particles (as for elec- 
trons or holes), that plays the central role in observables 
like current and noise. 

In the following, we will first analyze the scattering 
at a junction between a ^MM and a normal metal lead, 
which turns out to be also pairwise in the properly-chosen 
Majorana basis. In other words, the scattering at the 
junction can be significantly simplified in terms of both 
its physical picture and its mathematical parametriza- 
tion. Next we will apply the simplified scattering pic- 
ture to compute current and noise in several examples 
of chiral MF interferometers, including prototypes like 
Fabry-Pero1ji^, Mach-Zehnder-iii^ and Hanbury Brown- 
Twiss two-particle interferometers^. The emphasis of 
these calculations is placed on demonstrating the way 
Majorana scattering amplitudes enter the expressions for 
current and noise, and interpreting the results in a con- 
sistent manner. Most of the examples used in this paper 
have been discussed in the literature, therefore we will 
not dwell on detailed discussions about the implications 
of specific results. Nevertheless our approach shows its 
elegance in dealing with chiral MF interferometers and 
in understanding them. 



II. JUNCTION BETWEEN A CHIRAL 
MAJORANA MODE AND A NORMAL LEAD 

To address the question how MFs can be involved in 
charge transport, it is necessary to analyze how they are 
contacted with a normal (metal) lead which is indispens- 
able for measurements. In the chiral MF interferometry, 



where the propagating pathways for MFs can be devised 
to serve various particular interests, junctions between 
xMMs and normal leads become a universal, as well as 
crucial, ingredient and deserve a dedicated investigation. 
In this section we will focus on such a contact problem, 
leaving the propagation of chiral MFs to be handled in 
the next section with specific examples of interferome- 
ters. 

To this end we consider a x^M tunneling-coupled 
to a normal lead which contains N transmission modes 
around the Fermi energy and is connected to an electron 
reservoir with chemical potential /i at its far end^. The 
junction is described by the scattering theory adapted to 
deal with interfaces between normal and superconduct- 
ing systems^. Namely, the N transmission modes in the 
normal lead are artificially doubled into their electron 
and hole copies, which are then treated independently 
(but constrained by the PHS; see below) with different 
occupation functions 
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h(E) = 1 - fe{-E) = 



g/3(E+5M) + 1 ' 



(1) 

(2) 



where /? — l/Zc^T with fc^ the Boltzmann constant and 
T the temperature, Sfi = fj, — fig with /x^ the chemical 
potential in the superconductor. Throughout this paper 
we assume 6fi is much smaller than the superconduct- 
ing gap, /is is fixed by grounding the superconductor, 
and energy E is always measured relative to fig- In or- 
der to remove the artificial doubling from actual physical 
effects, it suffices to count only contributions by states 
with E>0. 

Specifically, the scattering matrix, S, for the junction 
at energy E is defined by 



'j(+){Ey 

^^e^\E) 



SiE) 



^pir\E) 



(3) 



where 7 is the annihilation operator for chiral MFs, 

'tpe/h = {'iple/h,llj2e/h, ■■■,1pNe/hV IS the iV-VCCtor of 

electron/hole annihilation operators in the normal lead 
and the subscripts (+/— ) always stand for the outgo- 
ing/incoming states. The operators here are each nor- 
malized by (the square root of) the velocity of the cor- 
responding mode, such that S" is a (2A'^-|- l)-dimensional 
unitary matrix as a result of quasi-particle probability 
current conservation. We emphasize that S by itself does 
not respect charge conservation because of the implicit 
presence of the TSC upon which the ^MM has to live. 

Often it turns out to be more convenient to work with 
the scattering matrix in an alternative basis -the basis 
is called the Majorana basis and hence the scattering 
matrix is denoted by S]\i- which is defined as follows 



(+)(T^\ -^*n^) I (-) 



r]'^+>{E) 



ri'-KE) 



(4) 
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where 



(5) 
(6) 



with Ijv the A^-dimensional identity matrix acting within 
the space of electron or hole modes. Evidently Sm is also 
unitary, and is related to S* by a change of basis in the 
normal lead: 



S(£) = 



The inherent PHS in the current problem imposes that, 
by definition, 



(8) 
(9) 



both of which apply to outgoing and incoming states sep- 
arately. The latter of the above equations implies 



(10) 



which allows us to regard the scattering modes repre- 
sented by //-operators as artificial Majorana modes. The 
scattering matrices constrained accordingly by the PHS 
satisfy: in the electron-hole basis. 



^ ^S*{E)(^ ^ \=S{-E), (11) 



where E^; = Ijv with the Pauli matrix inter- 
changes the fiavors of electron and hole; in the Majorana 
basis. 



Sl,{E) ^ Sm{-E). 



(12) 



It is clear that at i? = 0, Sm becomes a real orthogonal 
matrixes. 

Indeed, in the rest of this paper, we consider the low- 
energy limit where the energy dependence of the scat- 
tering matrices in Eqs. ([TT|) and ([T^ can be dropped 
(accordingly we will also drop the energy dependence of 
the operators for the scattering modes). As a result Sm 
is always taken to be a real orthogonal matrix with de- 
terminant 1 (i.e. Sm G S0{2N + 1)). This significantly 
simplifies the theory that will be developed in the follow- 
ing. 



A. General strategy 

In this section our goal is to understand the contact 
between a xMM and a normal lead by decomposing the 
scattering matrix at the junction into its physical part 
and its (physically) irrelevant part. The physical part is 



what is relevant to physical observables in the scattering 
process, and can be taken as the canonical form of the 
scattering matrix; the irrelevant part by definition does 
not make physical contributions. To achieve this goal 
we will work exclusively with scattering matrices in the 
Majorana basis, 5*^, in the current section. 

The advantage of working in the Majorana basis is that 
Sm^ being a real orthogonal matrix in the low-energy 
limit, is both a faithful and a convenient representation 
of the scattering processes which involve MPs. Such a 
representation naturally includes all possible scattering 
events at the junction without additional constraints as 
Eq. (fTTj) (i.e. the set of 5a/ in the low-energy limit is 
precisely SO{2N + 1), while the corresponding set of S 
is only a subset of SU{2N -t- 1)). This greatly facilitates 
the decomposition that we will perform in this section. 
The disadvantage of working in the Majorana basis is 
that the "artificial" Majorana modes in the normal lead, 
defined in Eq. ([5]), are not associated directly with the 
distribution of electrons/holes in the reservoir. Therefore 
we need to return to the electron-hole basis in the next 
section to compute physical observables such as average 
current and noise. Nevertheless we will see that the de- 
composition of Sm pays off not only in simplifying the 
computations, but also in understanding physically the 
results. 

The idea to decompose a scattering matrix into its 
physical and irrelevant parts is to notice that the physics 
of a single scattering event is unchanged under U{N)- 
transformations of the (electron) basis for the N degen- 
erate modes (separately for outgoing and incoming ones) 
in the normal lead. Decompositions following this idea 
are conventionally referred to as polar decompositions of 
a scattering matrix when the physical part contains only 
the transmission eigenvalues^^i^. In the current prob- 
lem the same idea gets slightly more complicated by the 
demand of the PHS, which insists that the valid trans- 
formations are of the form 



(13) 
(14) 



with V,W e U{N). In the Majorana basis, if Eq. (O is 
also imposed on the transformed operators (with tildes) , 
this means 




(+) 



where 



V^Un 
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(-) 



(15) 



W*I^N (16) 
(V,W e SO{2N)). 
Accordingly the decomposition of is given by 
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FIG. 2: Illustrations of the scattering at a junction between a 
single-mode (A'^ = 1) normal lead and a xMM, in the electron- 
hole basis of the normal lead (a), and in the properly-chosen 
Majorana basis (b). In the latter basis, only one Majorana 
mode in the normal lead is coupled to the xMM, the other is 
simply reflected with amplitude 1. 



When V and W are properly chosen, Sm contains a mini- 
mal number of parameters^ that are relevant to physical 
observables, and hence is the canonical form of Sm- 

As a consequence of the PHS, the decompositions here 
do not in general bring Sm into a diagonal form as con- 
ventional polar decompositions do (the diagonal form is 
even impossible in the electron-hole basis because the 
XMM cannot be coupled only to a single electron/hole 
mode at all). Instead the decompositions of Sm turn out 
to take the form of Euler decompositions: mathemati- 
cally this means Sm appears to be a combination of un- 
coupled (2D-)planar rotations; physically this means the 
scattering of many Majorana modes is essentially pair- 
wise. In the rest of this section, we will derive such 
an Euler decomposition first in the simplest case when 
iV = 1, and then in the arbitrary A'^-channel case. 



B. The iV = 1 case 

When TV = 1, Sm is of dimension 3. Namely, Sm 
is equivalent to a rotation in three dimensions. Such a 
rotation can be parametrized by three parameters, known 
as the Euler angles, as follows 



Sm — 
where 



1 W R[d) 
R{a) j I 1 



1 




(18) 



i?(e)^ f'^'J 

\ sm4 cos^ 



(e = a,/3,0) (19) 



is a planar rotation by an angle ^. 

The above decomposition is exactly what is dictated 
by Eq. (ITTl) . To see this, simply notice that 



(20) 



thus we can identify V = R{—a) and W = R{I3). Namely, 
if we redefine the electron and hole operators by trivial 




FIG. 3: Illustrations of the scattering at a junction between a 
multi-mode (A^ > 1) normal lead and a xMM, in the electron- 
hole basis for the normal lead (a), and in the properly-chosen 
Majorana basis (b). In the latter basis, only one Majorana 
mode in the normal lead is coupled to the xMM, the others 
are reflected in pairs and one is reflected with amplitude 1. 



gauge transformations (in the normal lead) 



(+) 




and redefine 17 according to Eq. ([S|), then the Majorana 
scattering matrix contains one single parameter and takes 
the form 



Sm-^ Q 1 j 



(22) 



Eq. ([22|) is the canonical form of the Majorana scat- 
tering matrix when iV = 1, which says that, effectively, 
the xMM 7 is coupled to one (artificial) Majorana mode 
fji in the normal lead, leaving the other (artificial) Majo- 
rana mode ?72 to be completely reflected with amplitude 
1 (see Fig. ^. 



C. The arbitrary A'^ case 

When > 1, there is no decomposition like Eq. ^TE\\ 
readily available. Nevertheless one can show that (see 
Appendix M . the (2A^ -I- l)-dimensional Majorana scat- 
tering matrix can always be decomposed as Eq. (|17|) 
with 



Sm — 



N 



(23) 



where R{9i) couples the modes 7 and 7)1, R{9i) with i > 1 
couples the modes fji and fji+N^i, and the last diagonal 
entry 1 signifies a complete refiection of mode fi2N ■ Note 
that by definition (see Eq. (O), fji and fji^jsf are partners 
derived from the same normal mode ^ie/Ui that is. 



m = 



1 
71 



(V'l. 



Vi+N 



(24) 
(25) 
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The concept of such pairs of Majorana modes turns out 
to be fundamental in relating Majorana scattering ma- 
trices to charge transport, as will be revealed in the next 
section. 

The above canonical form of Majorana scattering ma- 
trices should be no surprise after we have understood the 
iV = 1 case: the xMM 7 effectively couples to a single 
Majorana mode 771 - determined by the specifics of the 
junction - in the normal lead, which leaves its partner 
?7i+Ar - determined solely by 771 - to be coupled to the 
rest of the Majorana modes in the normal lead; fji+N 
then "picks" another (single) Majorana mode 772 to cou- 
ple to, leaving in turn 772+Af behind; this process goes 
on until the last one 772^ is left completely alone (see 
Fig. [3]). Interestingly, this re-pairing picture is a good 
analogy to the re-pairing of MFs that happens at the 
topological phase transition in the Kitaev model for ID 
p-wave superconductors^ - here, the re-pairing of Majo- 
rana modes happens when a normal lead makes contact 
to a 2D p-wave superconductor with a chiral Majorana 
edge mode. 

III. CURRENT AND NOISE FOR SPECIFIC 
EXAMPLES 

In this section we compute the average current and the 
zero-frequency noise power in various prototypical chiral 
MP interferometers. 

To this end we introduce for convenience two matri- 
ces: F, which is a diagonal matrix containing the Fermi 
distribution functions for all incoming modes, and (^2)^, 
which is another diagonal matrix containing the weights 
for charge current carried by electron and hole modes in 
lead ly (1 for electron modes and —1 for hole modes; 
for modes that do not belong to lead v). Explicitly, we 
define 

K 

^ = 0^^., (26) 
i/=i 

(E2), = (0 0,.) ® (a, (g) IatJ, (28) 

where K is the total number of contacts, f,^e/h is the 
Fermi distribution function for electrons/holes in con- 
tact v, is the number of transmission modes (without 
double counting) in lead v, 0^' is the empty matrix cor- 
responding to the modes in lead v' . 

We will use the following formulas from the scatter- 
ing theory adapted for normal-superconducting hybrid 
systems^Sii^: the time-averaged current at normal con- 
tact a, {la), is given by 

la^l I dETY[FAa\ ; (29) 
Je>o 



the average current flowing through the grounded super- 
conducting contact is given, owing to current conserva- 
tion, by (omitting the minus sign) 

Is^ J2 (30) 

QGnormal 

the zero-frequency noise power, defined as Pafj ~ 
J^^dt^{SL{t)Sii3{0) + Sif3{0)Sia{t)) with the current 
fiuctuation 6Iu=a,/3(t) ~ Iu{t) — (lu), is given by 

-Pa/3 = ^ / dETT[FAa{l - F)Ap] . (31) 
" Je>o 

In the above formulas we have used a shorthand notation 
defined as follows 

(S,),-5t(I],),5. (32) 

Note that S here is the full scattering matrix in the 
electron-hole basis, which includes K contacts but does 
not exphcitly include the xMM (see, e.g., Eq. (^5]) in 
Sec. nil B[) . The full Majorana scattering matrix Sm is 
related to S by 

K 

S = uISmU^, C/e = 0«7jv.. (33) 
1^=1 

Before we proceed presenting the results, we confirm 
that the transformations in Eqs. (IT3l) and (|14p indeed 
lead to no change in physical observables. That is, we 
verify that and P^^ are unchanged upon replacing S 
by S, given by 

S^V^SW^, (34) 
Ve^® y.] , W^e ^ e f . (35) 

The proof is straightforward after noticing that 

[W(s,F] = and Vi. : [V^, (S,),] = 0, (36) 

where [, ] stands for the commutator. Therefore we will 
adopt only the canonical form Sm in the actual calcula- 
tions, and also suppress the tildes from now on. 

A. A single normal lead probing a chiral Majorana 
loop: N = 1 

We start with the simplest setup: a single normal lead 
with a single mode is probing a closed xMM, or, a chiral 
Majorana loop; the interior of the TSC - enclosed by the 
chiral Majorana loop - may accommodate riy vortices, 
each allowing for one MBS at zero energy. This setup, as 
well as the setup that will be discussed in the next sub- 
section (Sec. IIIIBI) . has been discussed by Law, Lee and 
Ngi^. Here we reformulate the problem in a consistent 
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where 



FIG. 4: Effective scattering for setups with a single-mode 
(A'^ = 1) normal lead probing a chiral Majorana loop: the left 
panel case (a) can be extended into the right panel one (b), 
to model, e.g., the presence of an additional point contact or 
scattering at two coupled MBSs (see text for details). 



Majorana language, with extensions that are intriguing 
in their own right. 

We assume that the MBSs are all deep inside the TSC, 
such that the coupling between the MBSs and the xMM 
or the normal lead, which exponentially decays with re- 
spect to their spatial separation, is negligible. With this 
assumption the vortices are taken into account only as 
an additional phase (f> — nyW that the chiral MF picks 
up after moving in a full circle, and the acquired total 
phase hence is if = UyT: + n + 2Tr E/Ei, where the single 
TT comes from the Berry phase, and Ei = hvAi/L, with 
L the circumference of the loop, vm the group velocity 
of the xMM, is the level spacing for the chiral Majorana 
loop. 

Following our discussion in Section III Bl the effective 
scattering process for the current example in the Majo- 
rana language is pictured in Fig. 2^, and the Majorana 
scattering matrix is given by 



Sm — 



ri 
r2 



(37) 



where r2 = 1, and 



1 - roe"^ 



(38) 



with ro the (local) reflection amplitude of rji at the junc- 
tion connecting to the still unclosed 7 (namely, cos^ in 
Eq. (|22p ). Note that the above Majorana scattering ma- 
trix is implicitly energy-dependent because of the energy- 
dependence of <p; it is unitary, but not necessarily real 
except at i? = 0. 

Inserting Sm into Eqs. p9| and ([3T|) . together with 
psp , we find the average current and the zero-frequency 
noise power (auto-correlator) to be 



I^ll dE[l-^{r^r*)\5f{E), 



E>0 



P ^ — 



dE 



E>0 



{2[i~n{nr;)]e{E) 

inrlfSfiEy}, 



(39) 



(40) 



6f{E) = UE) - h{E) 



sinh 



(41) 



cosh (3E + cosh ' 
e{E) = ,UE)[1 ~ ME)] + h{E)[l ~ ME)] 

= -kBT^[UE) + h{E)]. (42) 

Note that Sf{E) is positive (negative) definite if (5/i > 
{6fi < 0), zero at equilibrium, while Q{E) is positive 
definite at finite temperature. As a convention, we shall 
intentionally keep the formal complex conjugate of real 
variables (e.g. r2) for a reason that will be clear in a 
moment. 



1. Current as interference of Majorana modes 

Let us now examine the expressions for the current Eq. 
(|39l) and noise Eq. (|40| closely. The net current in Eq. 
(|39l) is evidently attributed to the Andreev process: an 
incoming electron or hole, weighted effectively by the dif- 
ference of the Fermi distribution functions Sf, is reflected 
with "probability" 3fi(rir2) G [—1,1], and the missing 
charges are absorbed by the superconducting condensate 
as Cooper pairs. Here, negative 3?(rir2) means that a 
particle with an opposite charge with respect to the in- 
coming particle is reflected with probability |3'i(fir2)|. 
The noise in Eq. (|40|) has two contributions: the flrst 
term in the curly brackets gives the thermal noise which 
persists at equilibrium at flnite temperature; the second 
term gives the shot noise which vanishes either at equi- 
librium, where /e — fh, or upon total normal or Andreev 
reflections, where ^(rira)^ = 1 — 3fJ(rir2)^ = 0. 

The remarkable expression for the reflection "proba- 
bility", 3?(rir2), highlights one essential feature of MF 
transport, namely, charge is always "carried" by the in- 
terference between (a pair of) Majorana modest, even 
though an individual Majorana mode is charge-neutral. 
This can be heuristically understood as follows. Take 
one (quasi-particle) scattering state given by, in the Ma- 
jorana basis, \^{E)) = -^[ri\T]i{E)) + ir2\m{E))]; in 
the electron-hole basis, the same state reads \'i'{E)) = 
ae\ipe{E)) + ah\iph{E)) with Oe = (ri -|- r2)/2 and an = 
(ri — r2)/2; we immediately see that the current den- 
sity carried by this state is I{E) — (e//i)(|aep — la^P) = 
(e//i)5R(rir2), which is precisely the interference term be- 
tween the two Majorana modes. When ri and r2 are 
identiflcd with the reflection amplitudes in Sm, 
becomes the reflection "probability" that appears in Eqs. 
([55)1 and (Uni)- In order to stress the central importance 
of interference between Majorana modes in this work, we 
will on most occasions keep formal expressions similar 
to 5R(rir2) without simplifying expressions written with 
formally complex conjugate but in reality real variables. 
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FIG. 5: Differential conductance, dl/dV in units of /h, as 
a function of ^ = 2neV /Ei + (riv + 1)tt for tlie setup in Fig. 
U^. The resonance peak for ip equal to an even multiple of it 
and the (anti-resonance) valley for ip equal to an odd multiple 
of TT is a robust feature regardless of the coupling strength, 
parametrized by ro, between the normal lead and the xMM. 
This plot is essentially a reproduction of the results of Law et 
al.iS. The purpose of this reproduction is to facilitate com- 
parisons with our later results. 



2. The zero-temperature low-bias limit 

At zero-temperature and low bias {Sfi — eV <^ Ei), 
~ (n„ -I- 1)tt, ri ~ ±1, and Eqs. and (gHl) can be 
reduced to^^: 



for even n„: / = P = ; (43) 
for odd n„: I=—V,P^O. (44) 



The Fano factor, defined as F — P/{eI), is give by 1-l-ri 
and takes the value 2 if riy is even, or if is odd. 

Physically the two situations here are realized exactly 
when the chiral MF self-interferes maximally destructive 
or constructive upon reflection, hence a stationary state 
is disallowed or allowed at zero energy, ri — ±1 therefore 
corresponds to the reflection of a Majorana mode that is 
completely off- resonance (ri = -|-1) or on-resonance (ri — 
— 1) with such a stationary stated. As a consequence, the 
scattering process is a total normal reflection (/ = 0) or 
a total Andreev reflection (/ = 2e'^V/h) - in either case, 
the shot noise vanishes. 

For more general cases beyond the low-bias limit, 
where if is allowed to vary with energy continuously be- 
tween an even and an odd multiple of tt, we plot the zero- 
temperature differential conductance dl/dV, as a func- 
tion oi ip — 2i:eV/ El + [n^ + l)7r and in units of e^//i, in 
Fig. [5] with different tq. It is clearly seen that a quan- 
tized resonance peak is established at (/? = (mod 27r) 
and completely suppressed at = tt (mod decreas- 
ing ro broadens the resonance peak but does not remove 
such a general feature. 




FIG. 6: Differential conductance, dl/dV in units of e'^ /h, as 
a function of ip = 2neV/Ei -\- n -\- tf) with (f) = n„7r for the 
setup in Fig. |4j3. We set = in the upper panel and 
set = (/>' — TT in the lower panel; for both panels L = L' 
is assumed and ro is fixed to be 0.7. In the upper panel the 
resonance peak splits into two with their distance proportional 
to the energy splitting of the two coupled chiral Majorana 
loops; in the lower panel the resonance peak persists whenever 
p is an integer multiple of n as long as the coupling between 
the two chiral Majorana loops is finite. 



3. Extended setups 

The present simple setup can be slightly extended by 
including another chiral Majorana loop, labeled by 7', 
which couples only to the original loop 7 with tunneling 
amplitude t (see Fig. |4}d). This extended setup can be 
used to model, for example, the effect of an additional 
point contact formed by confining the path of the orig- 
inal xMM, or scattering upon two coupled MBSs. The 
computation for the extended setup requires a minimal 
effort, which amounts to replacing ip in Eq. p8|) with 
(/?oxt, defined as 



r — e^'f 

<^cxt + arg(- — ), 

1 — re^f^ 



(45) 



where <p' — EL' /Hvm + tt -\- (p' with <p' the magnetic flux 
inside 7', and r ~ y/l — is the local reflection ampli- 
tude between the two chiral Majorana loops. The ex- 
pressions for the average current and the noise power are 
unchanged from Eqs. dSQ]) and pO)) . 

For simplicity we consider (p — (p' — tt and L = L'. 
When the coupling between the two loops are turned off, 
we know from previous discussions that a resonance-led 
total Andreev reflection occurs at _B = (more generally, 
at E = nEi with n an integer; we will focus on the energy 
range around zero) . Upon turning on the coupling - even 
with infinitesimal t - the resonance at zero energy is im- 
mediately switched off because t^cxt has a tt shift from p 
and hence ri changes sign. Indeed, since the energy spec- 
trum of the coupled loops is given by e^'^^xt = 1 (which is 
also equivalent to the resonance condition ri = ~1), we 
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FIG. 7: Majorana scattering channels for Fabry-Perot inter- 
ferometers considered in Sec. IIII Bl the upper panel case (a) 
can be extended to the lower panel case (b) to model the pres- 
ence of an additional point contact or two coupled MBSs (see 
text for details). 
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find that the resonance shifts to i? ~ ±(i?;/27r)t when 
t ^ 1. This is in close analogy to a system in which two 
MBSs (e.g. at the two ends of a chiral p-wave supercon- 
ducting wire) are coupled and lifted from zero energji^l. 
It deserves to be mentioned that if the current case is 
modified by setting one of the phases 4> or (p' to be and 
the other tt, then the resonance occurs at zero energy 
regardless of the coupling strength between the loops as 
long as it is finite. Differential conductances for both 
situations, — = or tt, are plotted in Fig. [Sj 



B. Fabry-Perot interferometer: = 1 

Our second example deals with the Fabry-Perot in- 
terferometer composed of two single-mode normal leads 
(labeled by L and R) probing a chiral Majorana loop si- 
multaneously. We assume that the two normal leads are 
mediated only by the xMM and other tunneling processes 
are negligible. Fig. [7^ illustrates the effective scattering 
in the Majorana language, and the corresponding Majo- 
rana scattering matrix is given by 



Sm = 



(TLi Ilb. 

rL2 

^RL Trx 

V TR^, 



where = ''i?2 ~ 1, and 



TRi = Z{rRQ - rioe*'^), 
tLR = —ZtLQtRae^'^ , 

tRL = —ZtLotRO, 

Z ={l- rLoTRoe' 



(46) 



(47) 
(48) 
(49) 
(50) 
(51) 



with the subscript meaning that the reflec- 
tion/transmission amplitude is local at the junction (cor- 
responding to Eq. (1^ . ro = cos6', to = sin 9, for two 



FIG. 8: Differential conductance measured at different con- 
tacts {L,R,S) and in units of e'^/h, as a function of <^ = 
2TTeV/Ei + (uv + l)7r for the Fabry-Perot interferometer in 
Fig. [7^, where vlo is fixed to be 0.7 and the two normal con- 
tacts are equally biased. In general, resonance peaks (anti- 
resonance valleys) of the differential conductances ai (fi — 
{tp = ±7r) persist as a robust feature. Asymmetric coupling 
at the two junctions {tlo tro) results in unbalanced reso- 
nance peaks in dlh/dV and diR/dV , while the resonance peak 
in dis /dV remains quantized at 2e^ /h. Particularly when one 
of the junctions is shut off (e.g. rRo = 1), the result returns 
to what is shown in Fig. [5] where the resonant- Andreev- 
reflection is completely local. This plot is again essentially a 
reproduction of the results in Ref.^'"*, presented here for the 
purpose of comparison with later results. 



leads separately) . Note that we have chosen a gauge such 
that ip belongs entirely to the lower arm. 

The average current and the zero-frequency noise 
power are given by 



II - 
Pll 



(52) 



_E>0 



/ dE[l-^{rLirl^)\5fL. 
Je>o 

dE{2[l~^{rLirl^)]eL 

+ \Tlr[Qr~Ql + \ ULa~fRb)^]], (53) 

Plr = Prl 

I dE{-]-)^{tLRrl^tRLT*i,^)5fL5fR, (54) 

where T^r = t^^ptLR- Ir and Prr can be obtained by 
interchanging subscripts L and R in the expressions for 
II and Pll. 

In the presence of multiple normal contacts, it is par- 
ticularly interesting to look at the average current Is 
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flowing through the grounded superconducting contact. 
Assuming the two normal contacts here are equally bi- 
ased {S^L ^ Sfij^ = eV), Is is simply given by 

Is^Il+ Ir 

(55) 

where we have used the explicit expressions of the reflec- 
tion amplitudes. An immediate observation based on the 
above result is that the resonance peak in the differen- 
tial conductance dIs/dV at (^9 = remains quantized at 
2e^//i regardless of the coupling of the two normal leads 
to the xMM (see Fig. H]). 



1. Cross-correlation as two-MF-exchange effect 



0.5 




FIG. 9: The kernel pLj?(i?) of the integral for the cross- 
correlator Plr in Eq. (|54|l . excluding the factor due to Fermi 
distribution functions, as a function of ip{E). The coupling 
between the normal leads and the xMM is assumed to be sym- 
metric. Situations with various coupling strength are plotted. 
Between resonance {ip = 0) and anti-resonance {ip = ±7r) 
points (cf. Fig. [8]), Vlr changes sign. 



Formally the expression (|52p for II (or Ir) is the same 
as Eq. ([5^ in the single lead case, meaning that one 
contact does not contribute to the average current at 
the other contact regardless of the electron/hole occu- 
pation in either contact. Naively this can be seen as a 
result of the charge neutrality of the xMM that bridges 
two contacts. A more precise interpretation, however, is 
that the Majorana scattering states sourced from differ- 
ent contacts (e.g. t/li transmitted from contact R and 
VL2) are not phase coherent when only single-particle 
scattering is considered, therefore cannot lead to finite 
average current according to our discussion in Sec. IIII Al 
This statement will be clear after we further examine the 
cross-correlation between the two contacts. 

The charge-current cross-correlator Plr (or Prl), 
given by Eq. (|54p. is manifestly non- vanishing at low 
temperature and low bias if, and only if, both contacts 
are biased {SfLM{E) 0) and exchange of Majorana 
states is allowed {tLRtRL 7^ 0). This is despite the fact 
that the mediating xMM carries no charge current by it- 
self. Indeed, the non- vanishing cross-correlation here is 
entirely due to those scattering events that involve two- 
particle exchange and permit coherence between scatter- 
ing states from different sources, and hence is entirely 
the "exchange noise". Such an exchange contribution 
is irrelevant for the average current but present in noise. 
Meanwhile, the absence of other types of noise (e.g. ther- 
mal noise, partition noise) in the cross-correlator again 
highlights the essential role that the interference between 
Majorana modes plays in charge transport. 

In contrast to the cross-correlation, the auto- 
correlation for contact L, given by Eq. (j53p . contains 
both the equilibrium thermal contribution (the first term 
in the kernel) and nonequilibrium contributions (the rest 
two terms). Contact R manifests itself in Pll only away 
from equilibrium (except for modifying r^o to r^i), which 
is clearly different from a conventional normal-metal- 
superconductor- normal-metal device^. 



2. The zero-temperature low-bias limit 

At zero-temperature and low bias {eV <^ Ei, equal 
for both normal contacts), the average current and the 
zero-frequency noise power can be reduced to^: 



for even 

1 e^V, 
lL = In = -^Is = —il-n), 

Pll = Prr = Plr = Prl = ^(1 - rf) , 



''LP + rRQ 
1 + rwrRo 



(56) 
(57) 
(58) 



ri = = tri 
for odd 71^: 

II = ^(1 ~ ri), Ir = ^(1 + ri). Is - —V, (59) 
n n n 

Pll = Prr = -Plr = -Prl = '^{l - rl) , (60) 

ri = TLi = -Tri = . (61) 

1 - rLorR.0 

The Fano factors, for individual normal leads defined 
as Fa = Paa/ i^Ia), and for the total noise power (or 
equivalently, for the superconducting contact) defined as 
P = (J2a pPc,0)/{eY.a^a), whcrc a,/3 L,R, are given 
by 

1 + ri 

for even n„: Fl = Fr ^ — - — , = n ; (62) 
for odd n,: Fl = Fr = F^O. (63) 

We see that the even/odd Uy cases represent rather dis- 
tinct processes. For even the conductance and Fano 
factors of both contacts L and R are equal. This indi- 
cates that the statistics of the electron transport is in fact 
identical at these leads and the process consists of a ran- 
dom sequence of crossed Andreev refiection processes^^. 
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FIG. 10: Differential conductance measured at the supercon- 
ducting contact, dIs/dV in units of e^/h, as a function of 
= 2iveV/Ei + TT + (f) with (f) = UvT^ for the extended Fabry- 
Perot interferometer in Fig. [TId. We assume (j) — (f)' , L — L' 
and r_Lo = J'flo = 0.7. Situations with various t between the 
two loops are plotted. When t = 0, this setup is simply two 
copies of the single-lead setup in Fig. |4l when t = 1, this setup 
reduces to the original Fabry-Perot interferometer in Fig. [7^, 
with ip + ifi' ~ n playing the role of ip in the original setup. 
For the case of = 0' — tt and L = L' , Is does not depend 
on t at symmetric coupling {rLo = '"flo), and its expression 
reduces to Eq. (|55p (see also Fig. [8j;) with ip substituted by 
2(/?. The plot shows that quite generally, dIs/dV peaks at 
resonance (upon eigenstates of the coupled loops) points, and 
is suppressed at anti-resonance points. 



In contrast for n.y odd, electrons are drawn at different 
rates from both contacts but such as to generate a per- 
fect stream of Cooper pairs entering the superconductor 
noiselessly {F — 0). Meanwhile, the cross-correlators 
change from being positive for even n„ to being nega- 
tive for odd , indicating increasing importance of local 
Andreev reflections. 

For a more general picture of the scattering processes 
beyond the low-bias limit, we plot in Fig. IH] the cross- 
correlation noise power spectrum (the kernel of the in- 
tegral in Eq. ([M]) excluding the factor due to Fermi 
distribution functions) Vlr = {~\)'^{tLR'''*]^2^RL'f*ji2) '"^ 
a function of ip. By using the explicit expressions of 
the scattering amplitudes, one finds clearly that Vlr 
changes sign between resonance (iy9 = mod 27r) and 
anti-resonance [ip — n mod 2tt) points (cf. Fig. |8]). The 
interpretation for this sign change is in fact the same as 
our preceding discussion for the low-bias cases. 
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FIG. 11: Cross-correlation noise power spectrum Vlr as a 
function of ip{E) for the extended Fabry-Perot interferometer 
in Fig. [Tja. In the upper panel we assume (j> = 4>' , and in 
the lower panel (fi = <j}' — n; for both panels we take L — L' 
and tlo = ruo = 0.7. Situations with various transmission 
amplitudes t between the two loops are plotted. Clearly be- 
tween resonance and anti-resonance points (cf Fig. llOp . Vlr 
changes sign. 



need to be made 

TLi = Z[rLo + rfioe^'^+^'^ - ^e'^^ + TLoVRoe''^')] , (64) 

TRi = Z[rm + rLoe'('^+^') - r(e*'^' + VLoTRoe''')] , (65) 

tLR = ZtLotmte'^'^+'^'\ (66) 

tRL = -ZtLotRot, (67) 



(68) 



((38|) at 
gZll-dnni) at r = 0, 



Note that the above amplitudes reduce to Eq. 
r = a/1 — = 1, and to Eqs. 
t — 1. We plot the differential conductance dis / dV and 
the cross-correlation noise power spectrum Vlr for the 
present case in Fig. [TU] and [TTl respectively. Similar to 
our previous discussions, the coupling between the two 
loops (assumed again to be of the same circumferences) 
leads to split resonance peaks if = 0', or a halved period 
a (f) ~ (p' — TT, for the differential conductance dIs/dV; 
the sign of Vlr shows consistently changes between res- 
onance and anti-resonance points upon eigenstates of the 
coupled loops. 



3. Extended setups 

As in the single lead case discussed in Sec. HIT Al the 
present setup can be extended in the same way by in- 
cluding another chiral Majorana loop which couples to 
the original one (see Fig. [7}d). For the extended setup, 
the following substitutions for the scattering amplitudes 



C. A single normal lead probing a chiral Majorana 
loop: N > 1 

Our previous two examples deal with single-mode 
{N — 1) leads, which are simple but nevertheless reveal 
most of the interesting physics. The real strength of the 
decompositions discussed in Sec. |lTl however, becomes 
obvious only in dealing with multi-mode {N > 1) leads. 
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FIG. 12: Effective scattering for a setup with a multi-mode 
(A^ > 1) normal lead probing a chiral Majorana loop. Only 
the (single) Majorana mode coupled to the chiral Majorana 
loop can "feel" the phase picked up in the loop. 



In this section we start again with a single-lead setup 
similar to that in Sec. IIII Al but with TV > 1 (see Fig. 

According to Sec. Ill C[ the Majorana scattering matrix 
for this setup can be written as 




FIG. 13: Effective scattering channels for a generic Fabry- 
Perot-type setup with K multi-mode normal leads coupled to 
a single chiral Majorana loop. In each lead, only one Majo- 
rana mode is coupled to the chiral Majorana loop, and hence is 
involved in transmission between contacts (any direct tunnel- 
ing between contacts is ignored); the other Majorana modes 
(omitted in this illustration) are locally, pairwise reflected. 



Sm — 



r2 -t2 
t2 r2 







V 







Tn -tN 

In tn 



(69) 



rN+iJ 



where r^+i = 1, ri is given by Eq. (|3S|) . and = cos6'.i 
and U = sin 9i with 1 < i < A'^ in accordance with Eq. 
([23l) . Note that the basis for this scattering matrix is 
reordered from the original one defined by Eq. ([5]), to 
be ?7 = {r]i,rii+N,V2,'n2+N,—,r]N,VN+N)'^, where one 
should keep in mind tliat r/j and ?7i+Ar are the partners 
associated with the i-th original mode as given by Eqs. 
dMl) and (123. 

In terms of the elements of this scattering matrix 
the average current and the zero-frequency noise power 
(auto-correlator) are: 



P = 



If di?^[l-5R(r,r:+i)]5/, (70) 

n- Je>o 

2 r ^ 

^ di?{^2[l-5i(r.r*+i)]e 
Je>o '^^^^ 

^ 1 

i=l 
N 

+ Y.mrr*-ikr*+i)Sf}. (71) 



does not appear directly in the average current similar to 
the case of scattering between contacts in a Fabry-Perot 
interferometer. The noise power, by contrast, contains 
not only the auto-correlation of each individual mode 
(the first two terms in Eq. ^T^i; cf. Eq. (gO])), but 
also the cross-correlation between different modes (the 
last term in Eq. cf. Eq. (|M)) ). 

In general, {r^ : 1 < i < N} are parameters deter- 
mined specifically by the details of the junction between 
the normal lead and the TSC, and do not necessarily lead 
to quantized conductance regardless of the value of if. 
This is certainly an important consequence of the pres- 
ence of multiple modes in the normal lead. An interesting 
exceptioni^, however, is the zero-temperature low-bias 
conductance when N — 2 and n„ is odd, where one finds 
ri = — r3 = —1, thus G — 2e^/h and P — 0. Wimmer et 
alii^ have proposed this exceptional case to be a robust 
transport signature for detecting a MBS, with which a 
zero-energy eigenstate of the chiral Majorana loop can 
be identified^. 



D. K multi-mode normal leads coupled to a chiral 
Majorana loop 



Clearly the total average current is a sum of contribu- 
tions from all A'^ modes individually: and r^+i are re- 
flection amplitudes for rji and rji^ m and their interference 
term 3?(rir*_|_i) registers the outgoing current contributed 
by mode i] scattering between different modes (i.e. ti) 



As the final example for the Fabry-Perot-type inter- 
ferometers making use of a single x^^i discuss a 
generic setup with K normal contacts each admitting Na 
(a = 1, K) transmission modes in the lead coupled to 
the xMM (see FigfO)). The Majorana scattering matrix 
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for this setup is given by 

/ Sll Sl2 • ■ ■ SiK \ 



Mach-Zehnder 



Sm = 



S21 S22 



S2K 



with 



(72) 



\SKl SK2 ■■■ SkkJ 



r2{a) -t2{a) Q 





where 

?-l(a) 





Ziroia) 



{a ^13) 

n ^o(/3)), 



^JV„ + 1(q)/ 

(73) 
(74) 



(1 



n^o(/3)) \ 



(75) 

(76) 
(77) 



and Ti^a) and ti(^a) with « > 1 take the values according 
to Eq. (|23p for contact a. The subscript again means 
that the reflection/transmission ampHtude is focal at the 
junction and contains no (/^-dependence. The basis here 
for each contact is ordered the same way as in (|69|) . 

The average current and the zero- frequency noise pow- 
ers are given by 



(78) 



p(0) 



dE- 



E>0 



(79) 



a.h—e.h 



J E>0 ^ 



(a ^ /3) (80) 



where P^q takes the form of the single-lead auto- 
correlator ((7T|) . 

A remarkable fact in these results is that in the pres- 
ence of many modes, the iy9-dependence of the average 
current or the auto-correlators, carried only by ''i(a) and 
ia^, can be overwhelmed by the large iy9-independent con- 
tributions from rj(c) and t^^Q,) (i > 1) at imperfect in- 
terfaces between normal leads and the TSC; the cross- 
correlators, however, preserve the (/p-dependent terms as 
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¥u, Vui, luilui 

Hanbury Brown-Twiss 





FIG. 14: Illustrations of the scattering channels in a Mach- 
Zehnder interferometer (a,b) and a Hanbury Brown-Twiss in- 
terferometer (c,d), in the electron- hole basis for the normal 
part (a,c), and in the properly-chosen Majorana basis (b,d). 
In the latter basis, the Majorana scattering matrix for a junc- 
tion between a normal channel and two xMMs is simply an 
identity matrix. In both setups the outgoing leads are con- 
nected to electron reservoirs that are grounded. 



their sole contributions, hence will optimally manifest the 
information embedded in the transmission of xMM. 

Without performing further quantitative analysis for 
this example, we simply point out two recurrent mes- 
sages in the above general results: the average current 
at one contact is only contributed by the interference be- 
tween pairs of Majorana modes in the corresponding lead, 
and the current cross-correlation between two contacts is 
purely due to coherent exchange of two MFs sourced from 
the two contacts. These two messages together recapitu- 
late the crucial role that interference between Majorana 
modes plays in charge transport. 



E. Mach-Zehnder and Hanbury Brown-Twiss 
interferometers 

So far we have been focusing on the junctions between 
normal leads and a single xMM, and the Fabry-Perot- 
type interferometers which comprise only such junctions. 
Alternatively another type of junctions, composed of two 
xMMs with opposite chirality and connected to a chiral 
normal channel, can be used to build interferometers of 
the Mach-Zehnder typ o^^i^^ or the Hanbury Brown-Twiss 
type^" (see Fig. 

The junctions between two xMMs and a chiral normal 
channel are reflectionless and can be described by two- 
times-two scattering matricesi^iiS,^ where the two compo- 
nents at the normal side correspond to electron and hole 
modes or two Majorana modes depending on the basis. 
Assuming again the scattering matrix at the junction is 
energy-independent at low enough energy, the Majorana 
scattering matrix belongs to 5*0(2) owing to the PHS. 
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It follows that by employing the same strategy as pre- 
sented in Sec. Ill A[ the canonical form of the scattering 
matrix, in the properly chosen Majorana basis, is simply 
an identity matrix. In other words, the junction effec- 
tively splits or joins, without mixing, two xMMs. This 
certainly simplifies significantly our discussions which we 
present below in terms of specific examples. 



contact R. As a direct consequence, the cross-correlator 
in Eq. (|88p involves only the thermal fluctuation in con- 
tact L, while all the other contributions (e.g. partition 
or exchange noise) are absent due to the one-way nature 
of the present setup. 

2. Hanbury Brown-Twiss interferometer 



1. Mach-Zehnder interferometer 

The Mach-Zehnder interferometer illustrated in Fig. 
[Hk has been investigated by Fu and Kaneii, and 
Akhmerov, Nilsson and Beenakker— . Here we simply 
reformulate the problem in the Majorana language we 
have been speaking so far, to echo some general features 
of Majorana transport that have been discussed in the 
previous part of this paper. 

The Majorana scattering matrix for the Mach-Zehnder 
interferometer (see Fig. [Hb ). defined as 



VRi 

is effectively given by 
Sm 



S 



M 



ti 

t2 



(81) 



(82) 



where we can choose a specific gauge such that ti — 
with if = [n^ + 1)tt + 2ttE/Ei, and ^2 = 1- 

The average current and zero-frequency noise power 
can be easily obtained: 



Il = t dE6fL, 
Je>o 



dE ^{ht*) SfL, 



E>0 



Is = ll dE[l-^{t^tl)\5fL. 
" Je>q 



Pll^^ dEQ 



E>0 



(83) 
(84) 
(85) 
(86) 
(87) 



" Je>o 

Plr ^Prl^~^ I dE ^(titl) Ql , 

" JE>0 



where 3?(ii<2) = costp and 3(tit2) = sirup when the ex- 
plicit expressions of ti and t2 are taken. 

Clearly the outgoing current in Eq. (|84l) appears to be 
the interference between two ^MMs - the transmission 
"probability" ^{tif^) has a physically identical interpre- 
tation as the reflection "probability" 3?(rir2) appearing 
in Eq. ((39)) . An obvious difference, however, between the 
present setup and a Fabry-Perot-type interferometer is 
that contact L here indeed provides coherent sources for 
the MF pairs -instead of single ones- that are drained by 



The Hanbury Brown-Twiss interferometer shown in 
Fig. [Hb has recently been investigated by Striibi, Belzig, 
Choi, and Bruder^". Here we reformulate the problem, 
again in the Majorana language. 

The Majorana scattering matrix for the Hanbury 
Brown-Twiss interferometer (see Fig. [Till), defined as 



(89) 



/Vui\ 




/VLl\ 


VU2 


— Sm 


1L2 


VDl 1 




\riD2/ 




\VR2/ 



is effectively given by 



S. 



M 



(tuL 























toR. 


V 








(90) 



In this case we can choose a specific gauge such that 



and tun — toR — toL — 1- 



The average current and zero-frequency noise power 
are given by 



lo^^ll dE6U, ia^L,R) 

JE>0 

lp=0, {f3 = U,D) 

r 

P»a.^^ dE e„ , {a = L, R) 

" JE>0 

Ppp-^T f dE[heL + eR) 

■ Je>o ^ 
+ 1 E {fi = U.D) (94) 



(91) 
(92) 
(93) 



a^b—e^h 



Pud - Pdu 

c2 



= ^ f dE{^lMtuLthRtDRt*DL)SfLSfR, (95) 

and all other cross-correlators are zero. Taking the ex- 
plicit expressions, we have ^ituLti/ntDiitJjj^) — cosip. 

As pointed out by Striibi et al^, a remarkable fact 
in the present setup is that although the average cur- 
rent at the two drains U and D vanishes identically, the 
cross-correlation between the two contacts is in general 
non-zero. This fact, in the Majorana language, can be in- 
terpreted as follows: the MF pair in either drain are nec- 
essarily from different sources, and are therefore phase 
incoherent for single-particle scattering which the aver- 
age current measures - the incoherence between MF pairs 
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implies vanishing average current; two-particle scattering 
allows coherence between MF pairs in both drains, and is 
included in the measurement of noise power - the coher- 
ent injection of two electrons/holes from the two sources 
leads to non-vanishing cross-correlation between current 
in the two drains. 

The scenario here bears close analogy to the N = 1 
Fabry-Perot interferometer discussed in Sec. IIIIBl where 
the average current contributed by a foreign contact is al- 
ways zero but the cross-correlation between two contacts 
is in general non-vanishing. The underlying physics in 
both scenarios, as it has been for essentially all our exam- 
ples, boils down to the central importance of interference 
in the charge transport with Majorana modes. 



IV. SUMMARY 

In this work we emphasized a scattering approach us- 
ing a global Majorana basis. Both the normal state leads 
and the system which intrinsically carries chiral Majo- 
rana states are described in the same basis. As a conse- 
quence, from a purely technical particle scattering point 
of view, Majorana scattering problems are rather similar 
to scattering problems in normal mesoscopic structures. 
However, Majorana modes are neutral and carry nei- 
ther charge nor (electric) current. In contrast to normal 
scattering problems the charge and current operators are 
not diagonal in the Majorana scattering states. There- 
fore, unlike in normal systems, current is not determined 
by transmission probabilities of Majorana fermions. We 
have shown that charge and current appear due to in- 
terference of pairs of Majorana modes. For averages of 
single particle quantities, like charge and current, to be 
non-vanishing, interfering Majorana modes must neces- 
sarily have originated in the same contact. In contrast, 
in two particle quantities, like the current noise, two- 
particle exchange permits interference of two Majorana 
particles even if they originate from different leads. The 
physics of neutral excitations is certainly an important 
field of future research. It can be expected that the the- 
ory of Majorana particles represents an instructive ex- 
ample of this development. 
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Appendix A: Decomposition of a 
(2A'^ + l)-dimensional orthogonal matrix 

We prove in this appendix that for an arbitrary matrix 
On € SO{2N + 1), one can find a decomposition 



O 



N 



such that 



1 




K 



R{9 



N 



1/ 



1 

w 



(Al) 



cos 6i — sin 6i 
sin 6j cos 9i 



{i^l,...,N) (A2) 



and V,yy e S0{2N) satisfy 



V2. -V2.-1A, ii = l,...,N) (A3) 
W2z = W2.-iA, ii = l,...,N) (A4) 

where V„ (Wn) stands for the n-th row of V (W), and 
the 2Nx2N matrix A is defined as 



A = Iat (g) 



1 

-1 



-1 











V 



1 
-1 0/ 



(A5) 



Before we start the proof, we point out that V (W) 
satisfying Eq. (|X3)) (Eq. (|X4)) ) is nothing but V (W) 
given in Eq. (jl6p written in a differently ordered basis 
(the purpose of such reordering is to make the central 
matrix at the right-hand-side of Eq. (jAll) taking a com- 
pact block-diagonal form). This can be shown by using 
explicitly 



^ ^ [ilN -itN 



(A6) 



hence we have, from Eq. ([T 



V 



where V and W are iV-dimensional unitary matrices (cf. 
Eqs. (IT3)) and (fT4|) '). It is then straightforward to see that 
by rearranging, for i — I, N , the i-th row and column 
to the {2i — l)-th and the {i + A^)-th row and column to 
the (2z)-th, the rearranged V and W satisfy Eqs. (jA3p 
and 

Now we prove our main claim by construction. 
First we notice that On can always be written as 



_ , cos 6*1 (— sin^i)^ 



(A8) 
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where both v and w are normahzed 27V-dimensional 
row vectors and Xq is a 2Nx2N matrix. We assume 
I cos 6*1 1 7^ 1 (namely, the upper-left corner element of 
the matrix O^v is not ±1), such that both v and w are 
well-defined up to an irrelevant sign. 

Then we construct two matrices V^^^, W'^) G SO{2N) 
by setting 



(1) 



(1) 



(A9) 
(AlO) 



and choosing V- and W,p^ with 2 < i < 2A^ so that the 



orthogonality and Eqs. (|A3I) and (IA4p are all satisfied 
(which can be easily shown to be always possible). We 
find 



1 



1 







/cos 9i ~ sin 01 
sin 6*1 


V 



0\ 



(All) 



where Oat-i €! SO{2N — 1). In other words, we have 



1 



v(i);^^io 



1 



T = 



R{0 



Obviously this procedure can be continued for On~i, 
and then for On -2 etc., until we arrive at Oq, which is 
1. Meanwhile we obtain Eq. (jAll) with 



V = 



vw 

l2(Af-l) 





v(2) ' ' 
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(1) 



(A14) 
(A15) 



It remains to show that V and W defined above satisfy 
Eqs. ([M]) and (|X4)) . To this end we observe A^ = -1, 
hence Eqs. (|A3p and (|A4p are equivalent to 



(A16) 



[V,A]==0, [W,A] = 0, 



where ATi = V^^^XqW^^^^ . 

Since all three factor matrices at the left-hand-side of 
the above equation (jAlip belong to SO{2N + \), so must 
do their product. It follows that 



cos 9i 
On-1 



(A12) 



where [, ] stands for the commutator. We also keep in 
mind that, by construction, all factor matrices at the 
right-hand-sides of Eqs. (jA14p and (|A15I) satisfy the 
above commutation relations. This immediately verifies 
the validity of our construction, and thus completes the 
proof. 
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In the literature there is another type of junctions consist- 
ing of one unidirectional "normal" mode and two Majorana 



modes with opposite chiralit y ' . The scattering theory 
for this type of junctions will be discussed in Sec. IIIIEI 
Ref. [3^ suggests to call an ensemble of such random ma- 
trices ORE (circular real ensemble) to distinguish it from 
the more familiar CUE (circular unitary ensemble). 
Following our strategy, the lower bound of the number 
of physical parameters is given by Dim[5'0(27V + 1)] — 
2 Dim[i7(A^)] = A'^, which turns out to be the actual num- 
ber of parameters in Sm- 

Of course, the real container of charges here is the su- 
perconducting condensate, the interference between Ma- 
jorana modes simply determines how normal leads draw 
electrons/holes from the container. 

Indeed, such robustness is of a more fundamental origin, 
namely the Beri degeneracy^, which we will discuss in 
details in our following work, using the Majorana picture 
similar to this paper. 



